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ABSTRACT 

■ We present the first results from a study designed to test whether, given high-quaUty spec- 

' trophotometry spanning the mid-ultraviolet-optical wavelength regime, it is possible to dis- 

o : tinguish the metal content and star-formation history of individual elliptical galaxies with 

sufficient accuracy to establish whether their formation history is linked to their detailed mor- 
phology and position on the Fundamental Plane. 

From a detailed analysis of ultraviolet-optical spectrophotometry of the 'cuspy' ellip- 
' tical galaxy NGC 3605 and the giant elliptical NGC 5018 we find that: 1) optical spectra 

, with A > 3500 A may not contain sufficient information to robustly uncover all the stellar 

' populations present in individual galaxies, even in such relatively passive objects as elliptical 

I galaxies, 2) the addition of the ultraviolet data approaching A — 2500 A holds the key to estab- 

lishing well-constrained star-formation histories for these galaxies, from which we can infer 
Q ■ a formation and evolution history which is consistent with their photometric properties, 3) 

I despite the superficial similarity of their spectra, the two galaxies have very different 'recent' 

C/3 i Star-formation histories - the smaller, cuspy elliptical NGC 3605 contains a high-metallicity 

' population of age ~ 1 Gyr, and has a position on the fundamental plane typical of the product 

of a low-redshift gas-rich merger (most likely at z 0.08), while the giant elliptical NGC 
, , 5018, with a sub-solar secondary population, appears to have gained its more recent stars via 

' mass transfer / accretion of gas from its spiral companion, 4) despite these differences in de- 

tailed history, more than 85% of the stellar mass of both galaxies is associated with an old 
(9-12 Gyr) stellar population of near-solar metallicity. 

This pilot study provides strong motivation for the construction and analysis of high- 
quality ultraviolet-optical spectra for a substantial sample of ellipticals spanning the funda- 
mental plane. 

Key words: cosmology: observations - galaxies: evolution - galaxies: formation. 



1 INTRODUCTION 

Despite recent advances in observational cosmology, the formation 
and evolutionary history of elliptical galaxies remains a controver- 
sial subje ct. Originally thou ght to be relatively simple, coeval sys- 
tems (e.g. lFaber et alj 19771) . the development of hierarchical mod- 
els in the late 1980s, coupled with the apparent failure of searches 
for luminous young ellipticals at high redshift, led many workers to 
the viewpoint that most ellipticals have acquired the bulk of the ir 
final mass at relatively low redshift jKauffmann & Charlotlll99§) . 
In fact, empirical evidence has been presented both for for- 
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mation o f the bulk of the s t ars in present-day e l lipticals at high 
redsh i ft jSender et alJ[T993: iDunlop et alj 1 19961 : [McCarthy etai] 
|2004 ICimatti et al.ll2004 iRenzinil |2006|) . and for the formation 
of ellipticals from the merging of spirals at lower redshift (e.g . 
Kauffman 1996; Baugh. Cole & FrenkI [l996l : ISaugh et"ai] Il99^: 
iKauffman & Chariot,. 1998iV A seemingly anti-hierarchical picture 
has emerged, with the oldest stellar popu l ations residing in the most 
massive galaxies jPe Lucia et al.ll2006l : IRenzinil l2006l and refer- 
ences therein). 

Recently, however, this long-running debate has taken an in- 
teresting turn, with evidence presented that there exist (at least) 
two distinct classes of elliptical galaxy, which differ in detailed 
morphology and occupy different areas of the fundamental plane 
jKormendv & Bended [19961 : 1 Faber et ai] |l997h . Specifically, some 
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ellipticals have flattened, or 'core' central surface brightness distri- 
butions, whereas others have power-law, or 'cuspy' central surface 
brightness distributions. The cuspy ellipticals generally have rela- 
tively s mall effectiv e radii, inter mediate mas s, significant rotational 
support jLauer et"al.) 1995; Fab er et aljl997 ), and show little or no 
radio or X-ray emission ( Bender et alJI 19891) . In contrast, the core 
ellipticals are generally more massive, have la rge effective radii , 
are pressure sup ported with very little rotation jLauer et"al]| 19951 : 
iFaber et al]|l997,) . and often show stro ng radio emission with high 
X-ray luminosities teeuing etaljl999l) . 

Recent N-body merger simulations indicate that equal-mass 
mergers, or mergers of bulge-dominated systems regardless of the 
mass ratio, lead to slowly rotating systems with flattened central 
surface brightness distributions and 'boxy' isophotes, while gas- 
rich unequal-mass mergers pro duce faster-rotating power- law el- 
lipticals, with 'disky' isophotes iKhochf ar & Burkerdl2005h . It has 
been suggested that the cor e galaxies had the b ulk of their stars in 
place a long time ago (e.g. iBender et alj|l993) , whereas the dissi- 
pative mergers prod u cing c uspy ellipticals occurred relatively re- 
cently. iGenzel et alj ( 1200 ih have claimed to find direct evidence 
of the latter process in action, with their near-infrared study of 
ULIRGs revealing that these violent gas-rich mergers are heading 
onto the region of the fundamental plane occupied by power-law 
ellipticals. This form s an interesting contrast with the discovery by 
lDunlopetal] ( l2003h that the hosts of quasars lie in the region of the 
fundamental plane occupied by giant core ellipticals. 

Such an apparently anti-hierarchical picture of elliptical 
galaxy formation (i.e. in which the most massive objects are the 
oldest) has received general support from several recent statistical 
studies (e.g. [Xa uffman et al. 2003; Cald well. Rose & Concannoi] 
l2003l ; lHeavenset aU2004' ). and various theoretical arguments have 
been advanced to explain how such a situation could arise within a 
CDM universe in which the growth of the un derlying dark matter 
distribution remains hierarchical. For example. lSilk & ReesHl99§) 
argue that outflows driven by the massive and luminous black 
hole in a young massive galaxy could lead to the ejectio n of gas 
and te rmination of star-formation activity. More recentlv. iBinnevI 
( l2003h argues that, in the presence of sufficiently dense hot at- 
mosphere, cold infalling gas tends to be ablated and absorbed by 
the hot atmosphere, preventing any late star-formation activity be- 
yond that produced by the initial infall of cold gas. It has also 
been argued that the flat cores in the most massive ellipticals may 
be the result of mergers between two gala xies, each of which al- 
ready contains a sup er-massive black-hole (Begelman et al. 1980'; 
iFerrarese & Merrittj i2000; Ravindranath, Ho & Filippenko 2002). 
In this picture, cores could form in the final stages of the merger 
when the black holes form a binary. Stars are then ejected from 
the central regions as this binary hardens. Recent N-body simu- 
lations support this idea, producing massive core elliptic als from 
the m erging of bulge-dominated galaxies (e.g. Khochfar & BurkertI 
I2OO5I and references therein). 

Despite this wealth of observational and theoretical activity, 
what is still lacking is a detailed and convincing determination of 
how elliptical galaxy star-formation histories vary as a function of 
position on the fundamental plane. If this could be properly estab- 
lished, it would be of enormous assistance in disentangling the rel- 
ative importance of the different mechanisms for elliptical galaxy 
formation currently under discussion. 

That this has not been achieved despite, for example, the avail- 
ability of optical spectra f or large samples of ell ipticals from the 
Sloan Digital Sky Survey dEisenstein et al.|[2003l) . may seem sur- 
prising but is due to two well known problems. First, there is the 



problem of age-metallicity degeneracy, with attempts to determine 
the ages of stellar populations in elliptical galaxies from optical 
data being largely foiled by the fact that colours and absorption 
lines (e.g. Mg) are degenerate to compensating changes in age and 
metal content ( e.g. IWorthevIl 1994 iTrageJI 19991) . Second, although 
progress has been made with the availabilit y of higher-resolution 
model spectra (e.g. lBruzual & Charlolll2003h . there is the problem 
that with limited spectral range, the data may not allow separation 
of different stellar age components. For example, while the Balmer 
lines are only moderately sensitive to changes in metallicity, even 
small amounts of recent star formation can strengthen the Balmer 
lines out of all pro portion to the actual mass involved in the burst 
( lTrageretall200ll) . 

In this paper we present the results of a pilot study designed 
to test whether, given high-quality expanded wavelength coverage 
extending into the ultraviolet, these problems can be overcome. 
Specifically we aimed to test whether, with accurate spectropho- 
tometry spanning A ~ 2500 —+ 8000 A, the age-metallicity de- 
generacy can be lifted, and the full star-formation history of the 
elliptical disentangled. Various authors have recognised the poten- 
tial of near-ultraviolet data to clarify the star-formation history of 
ellipticals. In particular, Fanelli et al. (1990) found that mid-UV 
colour and absorption line indices have dependences on temper- 
ature and metallicity which are strong and di stinct from those de- 
rived l ongward of 3500 A dLotz et al J 120001 , see also). More re- 
cently, iDorman. O'Connell. & Roodd2003l) have carried out an ex- 
tensive exploration of the usefulness of the mid-UV spectral region 
for determining ages and abundances in old stellar populations, and 
emphasise the importance of working with low-resolution spectra 
because of the present inadequacy of high-resolution models. 

We sought, therefore, to construct well-calibrated ultraviolet- 
optical spectra of moderate resolution (i.e. 5-10 A) for elliptical 
galaxies by combining publically available HST ultraviolet spec- 
troscopy data with existing optical spectrophotometry. NGC 3605 
and NGC 5018 transpire to be good representatives of the cuspy 
and giant core classes of ellipticals respectively. 

The structure of the paper is as follows. In the next section we 
briefly summarise the known properties of NGC 3605 and NGC 
5018, and describe the available ultraviolet and optical spectro- 
scopic data which can be combined to produce a properly flux 
calibrated ultraviolet-optical SED suitable for modelling. Then, in 
section 3, we describe the model fitting that we have undertaken, 
exploring models of varying sophistication to determine how com- 
plex the star-formation histories of these objects could be, and what 
results can be regarded as robust. In section 4, we present the re- 
sults of the model-fitting process, and summarise what can be de- 
duced with confidence from this study about the star formation his- 
tories of these two elliptical galaxies. Finally, in Section 5, we con- 
sider the implications for elliptical galaxy formation, and discuss 
the prospects for what could be achieved with a more extensive 
study, before summarising our conclusions in Section 6. 

Throughout this paper we assume a flat cosmology with Ha = 
70kms"^Mpc"\ r2„ = 0.3, and f^A = 0.7. The age of the Uni- 
verse in this model is 13.5 Gyr. 



2 DATA 

NGC 3605 is an elliptical galaxy with redshift z = 0.002228, 
and absolute magnitude Mb = —18.48. It is a member of the 
galaxy group NGC 3607. The UV spectrum (2260 - 3250, and 
3250 - 4760 A) was observed using the Faint Object Spectrograph 



(FOS) on the Hubble Space Telescope jPonder et al]|l998h . The 
optical spectram (3520 — 7340 A) was observed with the FAST 
spectrograph at the F. L. Wh ipple Observatory's 1.5m Tillinghast 
telescope (Ja nsen et alj|200cil) . The FOS observation used a circu- 
lar l"aperture, and the optical spectrum was observed using a 3'x 
3"slit, scanned over a distance equal to half the blue minor-axis 
optical diameter 

NGC 50 1 8 is a giant elliptical with redshif t 2 = 0. 009320 and 
absolute magnitude Mb = -22.22 (Schweizer & Seitzer 1992.) . It 
is a companion to the spir al galaxy N GC 5022. Both an optical 
dSchweizer & Seitzej[l992l) and an HI jlCim et al.|[l988l) bridge to 
this galaxy have been observed. The UV spectrum for this galaxy 
(2260-47 00 A) was also o bserved using FOS on the Hubble Space 
Telescope jPonder et ai]T998,) . The optical spectrum (3709 - 8000 
A) was observed using the Boiler and Chivens spectrog raph at the 
ESO 1.52m telescope at La Silla teica & AIloi3ll987l) . The final 
combined spectrum has a total wavelength range 2260 — 8000 A. 
Again, the FOS observation used a circular l"aperture. The optical 
spectrum was observed with a 5"x 8"slit. 

Clearly, the apertures for the UV and optical sections are not 
matched in either case. We compensate for this in part by allow- 
ing the relative normalisation for each spectral segment to float 
independently where possible (see §3). If there has been a sig- 
nificant episode of star formation in the last few gigayears, one 
would expect the starburst to have occurred in the central regions, 
so that the smaller FOS aperture may contain flux from a higher 
proportion of young stars than the larger optical apertures. How- 
ever, as the UV region of the SED is likely to be dominated by 
the flux from young stars anyway (thermally-pulsing asymptotic 
giant branch st ars are not significant contributers to the flux at these 
wavelengths, ( Maraston 2005)), in this case, the leverage for the ra- 
tio of young-to-old stars comes mostly from the optical spectrum, 
which mitigates the effect of un-matched apertures. 

The star-formation rate, estimated from [OIIJ3727 emission, 
an error on the FOS fluxes of 6%, and the relationship derived by 
iGallagher, Hunter. & Bushousd j 19891) is ~ (1-a few) x 10"" Mq 
yr~^ for both galaxies, which is as one might expect in early-type 
galaxies such as these. 

The errors on the spectra were estimated for each wavelength 
bin from the standard deviation of the nearest nine points to that 
bin, where the bin size was the full width half maximum of the 
spectral observation (8 A for the UV spectra, 6 A for the opti- 
cal spectrum of NGC 3605, and 1 1 A for the optical spectrum of 
NGC 5018). After the estimation of the errors, the spectra were 
re-binned, where necessary, to either the resolution of the lowest- 
resolution section of the stellar population evolutionary synthesis 
model spectra (10 A), or, where the resolution of the data was > 
10 A, the models were binned up to match the data (with errors 
propagated in quadrature). 

The resulting spectral energy distributions for these two galax- 
ies are shown in Figure[T] 



3 MODEL FITTING 

We have u ndertaken our analy s is usi ng the most recently updated 
models of iBruzual & CharlotI ( l2003il . hereafter BC03. The age 
range covered by this model set is 0.01 - 14 GyiQ. The models as- 
sume that solar metallicity ( Zq) is 0.02, and the available metallic- 

^ The resolution of the age grid of the modelled stellar populations is such 
that the 14 Gyr population is the closest to the age of the universe. For all 
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Figure 1. The combined UV-optical spectra of the two elliptical galaxies 
NGC 3605 and NGC 5018. For NGC 3605, the spectra have been rebinned 
to a resolution of 10 A. For NGC 5018, the UV portion is at 10 A resolution, 
and the optical part at 1 1 A. The normalisation is arbitrarily adjusted for 
clarity. 

ities are: 0.02, 0.2, 0.4, 1.0, 2.5 7jq. The wavelength resolution of 
our data is higher than the outputs generated by the other modelling 
codes. Although our data could be rebinned to their resolution, our 
preferred method is to rebin the models and leave the data its least 
manipulated state. 

The galaxy spectra were fitted with synthetic stellar popu- 
lation models in a number of ways. I nitially, we have used the 
MOPED algorithm jHeavens et alJuOOQ) . which allows the initial 
assumptions about the star formation history to be unconstrained, 
but does not necessarily return the simplest model consistent with 
the data. We have then fitted each set of galaxy spectra with a sin- 
gle stellar population, for each available metallicity in turn, with 
age and relative normalisation allowed to vary freely. Finally, we 
have fitted a two-component model, in which the age, metallicity 
and relative mass fraction of the two components are allowed to 
vary freely. In the latter two cases, the whole parameter space was 
searched, with the best fit determined via minimisation. The 
entire available wavelength range longward of 2500 A for each ob- 
served galaxy (i.e. 2500 - 7340 A for NGC 3605, and 2500 - 
8000 A for NGC 5018) was included in the fits. The sky lines at ~ 
6600 A and ~ 7560 A were excluded from the fitting in all cases. 

3.1 Fitting witli the MOPED algoritlim 

In order to determine the optimal multi-population fit for each spec- 
trum it is necessary to allow as free a star formation history as pos- 
sible. A parameterization of 1 1 single stellar populations, spaced 
logarithmically in time, each with metallicity varying freely, was 
all owed, together with a one-parameter slab extinction based on 
the iGordon et alj ( l2003h LMC curve. We chose to use the LMC 
extinction c urve as it is current ly the most reliably c haracterised 
dust s creen jGordon et al.ll2003l) . Extensive testing in IPanter et al.l 
(2006!) has shown that for non star-burst galaxies the choice of ex- 
tinction curve does not affect the recovered star formation history 
at optical wavelengths. 

In order to fully assess the resulting 23-dim ensional likeli- 
hood surface it was necessary to use the MOPED FI lHeavens et al.l 

practical purposes the spectrum at this age can be taken to be equivalent to 
that of a 13.7 Gyr population. 

^ Massively Optimized Parameter Estimation and Data compression 



4 L.A. Nolan et ah. 



bOOd ; iReichardt et alj 1200 ih algo rithm coupled wi th a Markov 
Chain Monte Carlo search method jPanter et al.l2003h developed to 
deal with the high resolution BC03 models (Panter et al., in prepa- 
ration). For each galaxy the best solutions from a series of 100 
randomly-seeded, conjugate-gradient searches were used to seed a 
Markov chain of 10*^ steps. The individual search traces were ex- 
amined by eye to confirm full exploration of the parameter space 
and estimate solution convergence. Due to the speed of the MOPED 
algorithm this analysis takes of the order three minutes per galaxy 
on a standard desktop machine, allowing a thorough investigation 
of a grid of different relative normalisations of the UV, visible and 
IR spectra, in order to reproduce the star-formation history as ac- 
curately as possible. 

For both galaxies the individual sections were scaled to opti- 
mally match in the overlap region and the final spectra normalise d 
to unity at 5500A then analysed with MOPED IPanter et a02006h . 



4 RESULTS 
4.1 NGC 3605 

4.1.1 MOPED result 

The spectrum of NGC 3605 was fitted with the MOPED algo- 
rithm, and Figure |2] shows the bestfit result of searching the 23- 
dimensional parameter space with optimal relative normalisation 
of the UV and optical sections. The results for NGC 3605 are dom- 
inated by two populations with central ages of 1.1 and 12 Gyr and 
metallicity Z = 2.32 Z© and Z = 1.00 Z0 respectively. The older 
of the two components accounts for 87% of the total stellar popula- 
tion by mass, and the second component accounts for the remaining 
13% of the stellar mass. The dust attenuation associated with the fit 
is E{B ~V)= 0.00251, consistent with the overall picture of an 
old elliptical with no recent star formation. The reduced of this 
fit is 1.774. 



3.2 Single-metallicity model fits 

Each available, single-metallicity, instantaneous starburst model 
was fitted individually to the spectra of the two galaxies. Each spec- 
tral segment (three for NGC 3605, and two for NGC 5018, see Fig- 
ure [T]l was normalised to a mean flux per unit wavelength of unity 
across its length. For each galaxy, the fit was performed simultane- 
ously on its set of spectra, with the model age varying freely. Unlike 
the MOPED fits, here we are able to allow the normalisation to vary 
independently for each section (to compensate for any deviations in 
the flux calibration). The best-fitting age for each metallicity was 
found via minimisation. 



4.1.2 Single-metallicity model fit 

The results of the single-metallicity fits to the spectra of NGC 
3605 are listed in Table [T] For each metallicity, the free param- 
eters are the age and the overall normalisation, which is allowed 
to vary freely for each section to compensate for any inconsis- 
tencies in the flux calibration of the spectral segments. Figure |3] 
shows the best-fitting single-population model spectrum, superim- 
posed on the observed spectra. The single stellar population model 
struggles to match the shape of the spectrum in the second segment 
(3300-4700 A), and many of the detailed features. This failure is 
borne out by an unacceptably large value of reduced namely 
x'=2.7. 



3.3 Two population fits 

Finally, a two-population model was constructed, where the age, 
metallicity and fractional contribution by stellar mass for each com- 
ponent were allowed to vary as free parameters, in order to deter- 
mine whether more than one significant episode of star formation 
has occurred. We have 



F2pop,\ = const{Xifzi,\,ti + Xjfz^,\,tj) 

Here, F2pop,\.t is the model flux per unit wavelength in the 
bin centred on wavelength A which is the sum of the two single- 
metallicity model fluxes (per unit wavelength in the bin centred 
on wavelength A), fzi,x,ti and fz ,\.t , which have metallicities 
and ages, Zi,ti Gyr, and Zj,tj Gyr, respectively. The fractional 
contributions of the two components, Xi and Xj, are allowed to 
vary freely, and again, normalisation to a mean flux of unity across 
the wavelength range of each spectral segment was carried out af- 
ter the addition of the two unnormalised fluxes. As for the single- 
component model fits, the two-component fitting was carried out 
simultaneously for each set of galaxy spectra, with the relative nor- 
malisation allowed to vary independently. 

A minimum-x^ fit was used to determine the best-fit values 
of Xi, Zi,ti, Xj , Zj and tj . Again, the whole parameter space was 
searched, with the best-fitting parameter values corresponding to 
the point on the agei-agej-Zi-Zj grid with the minimum calculated 



4.1.3 Two population fit 

The results for the two-population model fits to the spectrum of 
NGC 3605 are presented in Table |2] and Figure |4] A young, high- 
metallicity (1 Gyr, 2.5 Z©) component is still found, but the fit is 
much improved by the addition of a dominant mass-fraction of old 
(10-14 Gyr), lower-metallicity (Z0) stars. Degrading the spectral 
resolution to 20 A or 40 A makes no difference to the recovered 
best-fit parameters, although xi decreases to 1.3 and 0.9 respec- 
tively. 

4.2 NGC 5018 

4.2.1 MOPED result 

The spectrum of NGC 5018 was fitted with the MOPED algorithm, 
and Figure |5] shows the bestfitting model resulting from searching 
the 23-dimensional parameter space. The results suggest that NGC 
5018 is dominated by a population in the age bin centred on 8.5 
Gyr, with metallicity Z — 1.50 Zq (94% by stellar mass), but has a 
small contribution (5% by stellar mass) from the age bin centred on 

2.3 Gyr, with a metallicity of Z = 0.79 Zq. There is a very minor 
(0.07% by stellar mass) population in the bin centred on 1.1 Gyr, 
with a solar metallicity. The dust attenuation associated with the fit 
is E{B — V) = 0.077, consistent with the overall picture of an old 
elliptical with no recent star formation. The reduced x^ of this fit 
is 1.39. 
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Figure 2. Best-fitting model (tliin line) to the optimally-spliced spectrum 
of NGC 3605 (thick line) predicted using MOPED (Heavens et al., 2000). 
The skyHnes have been masked out of the fit. The dominant population, has 
age =12 Gyr, 1 .00 Zq and a relative stellar mass fraction of 0.87. There is 
also a significant (13% by stellar mass) contribution to the flux from a 1.1 
Gyr, 2.32 Zq population. See section 14.1.1 I for discussion. 




3000 



4000 5000 
restframe A / 



6000 



7000 




Figure 3. Top: the best-fitting single-metallicity model (1 Gyr, Z = Zq, 
thin black line) superimposed over the spectra of NGC 3605 (thick black 
line). Bottom: the contours enclose the regions where Ax^< 18.4. See 
section 14.1.2 I for discussion. 
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Figure 4. Top: the best-fitting two-component model (thin black line) su- 
perimposed over the spectral sections of NGC 3605 (thick black line). The 
two component populations are also shown; the dominant population (Z = 
Zq, age = 14 Gyr, M/Mg^; = 0.84) is in red, and the secondary stellar 
population (Z = 2.5 Zq, age = 1 Gyr, M/Mg^jj = 0.16) is in green. Bot- 
tom: the contours enclose the regions where Ax^< 16.8. See section l4. 1.31 
for discussion. 



Table 1. Table of results for the single stellar population fits to the spectra 
of NGC 3605 and NGC 5018. See sections lTOl and liiTIl for discussion. 
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4.2.2 Single-metallicity model fit 

The results of the single-population fits to the spectrum of NGC 
5018 are presented in Table[T] Figure |6] shows the observed spec- 
trum, with the best-fit model spectra superimposed. The minimum 
reduced is y^= 3.0, from which we may infer NGC 5018 can- 
not be described by a single, simple stellar population. The best-fit 
model has an age of 11 Gyr, and solar metallicity. This does not 
allow the UV flux to be well reproduced. 
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Table 2. Table of results for the two-component stellar population fits to 
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Figure 5. The best-fitting model (thin line) to the optimall y-spliced spec- 
trum of NGC 5018 (thick hne) predicted using MOPED iHeavens et"ai] 
I2OOOI) . The skylines have been masked out of the fit. The dominant pop- 
ulation, has age = 8.5 Gyr, 1.5 Zq and a relative stellar mass fraction of 
0.94. There is also a significant (5% by stellar mass) contribution to the flux 
from a 2.3 Gyr, 0.79 population. See section l4.2. 1 I for discussion. 



4.2.3 Two population fit 

The results of the two-population fit to the spectrum of NGC 5018 
are presented in Table [2] and Figure |7] The addition of a second, 
younger component (3 Gyr, 0.4 Zq) significantly improves the fit, 
particularly shortwards of the 4000 A break, drops from 3.0 
for the single stellar population, to 2.1 with the addition of two 
extra parameters. Degrading the spectral resolution to 20 A or 40 A 
makes no significant difference to the recovered best-fit parameters. 

decreases to 1.4 and 1.0 respectively, and at 40 A resolution the 
fitting recovers a dominant mature population of 11 Gyr, rather than 
12 Gyr. 
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Figure 6. Top: The best-fitting single-metallicity model (thin black line; 
age =11 Gyr, Z = Zq) super-imposed over the spectral sections of NGC 
5018 (thick black line). Bottom: the contours enclose the regions where 
Ax^< 18.4. See section 14.2.2 I f or discussion. 

5 DISCUSSION 

5.1 Age and metallicity determination 

In this work, model fits to long-baseline spectra have enabled the 
age-metallicity degeneracy to be lifted. Comparison of Figures 4 
and 8 (for NGC 3605) and of Figures 7 and 9 (for NGC 5018) 
demonstrates the improvement possible in the simultaneous deter- 
mination of both age and metallicity when a wide wavelength range 
is used. In Figures 8 and 9 only the optical (longwards of 3800 A) 
spectra have been fitted. Without the leverage from the UV spectra 
it can be seen that the age of the dominant component in NGC 3605 
is almost completely unconstrained, although the young component 
is still detected. For NGC 5018, although a mature, metal-rich pop- 
ulation is recovered, without the extra information provided by the 
shorter-wavelength data (which are dominated by younger stars) 
the fitting process instead tries to compensate in matching the shape 
of the spectrum by choosing a metal-poor mature population. This 
clearly demonstrates the power of the long-baseline fit in determin- 
ing age and metallicity, both in constraining the ages and metallic- 
ities of the components, and in selecting the correct components. 

For NGC 3605, all the detailed fitting (BC03: single stel- 
lar population, two-component long-baseline and optical-only fit- 
ting) agrees that there is a young (< 3 Gyr), metal-rich compo- 
nent present, but only the long-baseline two-component fitting can 
constrain the age of the dominant mature population. The results 
from the full two-component fit are in excellent agreement with the 
MOPED result. 

In Table [3] we compare our one- and two-component fitting 
using the BC03 models with the best-fitting parameters recovered 
using models from different authors, namely those of iMarastonI 
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Figure 7. Top: the best-fitting two-component model (tliin black line) su- 
perimposed over the spectral sections of NGC 5018 (thick black line). The 
two component populations are also shown; the dominant population (age 
= 12 Gyr, Z = 0.4 Zq, M/M gal = 0.88) is in red, and the lesser popu- 
lation (age = 3 Gyr, Z = 0.4 Zq, M/Mgai = 0.12) is in green. Bottom: 
the contours enclose the regions where Ax^< 16.8. See section l4.2.3l t"or 
discussion. 
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Figure 8. Top: the best two-component model (thin black line) fit to the 
optical-only spectral sections of NGC 3605 (thick black line). The two com- 
ponent populations are also shown; the dominant population (age = 5 Gyr, 
Z = Zq, M/M gal = 0.88) is in red, and the lesser population (age = 
1 Gyr, Z = 2.5 Zq, M/Mg^j = 0.12) is in green. Bottom: the contours 
enclose the regions where Ax^< 16.8. See section IsTl for discussion. 



(l2005l) and the PEGASE models dFioc & Rocca-Volmerangel 
[l993). These models have a wavelength resolution of 20 A (lower 
than the BC03 models) and, although the age grid is the same for 
all three model sets, the available metallicities are Maraston ( 2005): 
0.02, 0.5, 1.0, 2.0 Zq, and PEGASE: 0.02, 0.2, 0.4, 1.0, 2.5 Zq. 
Given the differences between the models, the agreement between 
the best-fitting parameters recovered for NGC 3605 from the var- 
ious models is remarkable. The minimum x^is significantly im- 
proved by the addition of a second component, and in all cases, the 
models find a young (< 1 Gyr) secondary component, and a dom- 
inant old (> 12 Gyr), solar metallicity component. The main dif- 
ference is that the older PEGASE models do not select the highest 
available metallicity for the young component, but a solar metal- 
licity instead. However, these results show the robustness of this 
method to both choice of model and wavelength resolution. The 
observed 2MASS J-Ks and J-H colours are reproduced by the best- 
fitting stellar population models to < ± 0.15. This is a reasonable 
match, as one might expect some discrepancy between the mod- 
elled and observed colours, as the observed colours are those for 
the entire galaxy, whereas the models represent only those stellar 
populations residing in the central regions of each galaxy. 



iTerlevich & Forbed (l2002h estimated the mean (luminosity- 
weighted) age of the central regi on of NGC 3605, by fitting the 
older stellar population models of IWorthev & Ottavianil ( Il997t) to 
H/3 and [MgPe] line indices. They find a mean age of 5.5 Gyr, and 
a mean metallicity of 1.67 Zq. The high metallicity is consistent 
with the high metallicity of both our component populations, and 
the mean age, which may be considered an upper limit to the age 
of the most recent burst star formation, is also consistent with the 
presence of the young population that we find. 

For NGC 5018, all the detailed BC03 model fits agree that this 
galaxy is dominated by a mature (> 10 Gyr) metal-rich (> Zq) 
stellar population. The power of the two-component long-baseline 
fit allows us, in addition, to disentangle the lesser, younger (3—6 
Gyr), metal-poor (Z = 0.4 Zq) population. 

Again, we compare our results from the BC03 models to those 
recovered from the Maraston and PEGASE models in Table [3] 
Again, the agreement between the various models is remarkable, 
demonstrating the robustness of this method. Adding a second com- 
ponent significantly improves the fit, and all the models find a dom- 
inant old (> 12 Gyr) stellar population, at the highest metallicity 
available. The BC03 and Maraston models agree that there is a sub- 
solar metallicity, young (2-3 Gyr) component, although again, the 
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Figure 9. Top: the best two-component model (thin black line) fit to the 
optical-only spectral sections of NGC 5018 (thick black line). The two com- 
ponent populations are also shown; the dominant population (age = 14 Gyr, 
Z = 2.5 Zq, M/M gal = 0.86) is in red, and the lesser population (age = 
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Table 3. Comparison of results for the various stellar population model fits 
to the spectra of NGC 3605 and NGC 5018. See section lsTl for discussion. 
The agreement between the different models is excellent. 
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older PEGASE models differ, by selecting an old metal-poor popu- 
lation to fit the blue component. The observed 2MASS colours for 
NGC 5018 are reasonably well-reproduced by all the models, to 
< ± 0.17. 

However, using an independent method (co mparing the ab- 
sorptio n line indices Ca II and H5/A4045), and the lWorthevI ( Il994) 
models. iLeonardi & WorthevI \2QQ&i find that a 2.8 Gyr, Zq popu- 
lation dominates the spectrum of NGC 5018 at 4000 A. They claim 
an old, metal-poor population cannot be c ontributing signifi cantly 
to the integrated blue light. Our BC03 and iMarastoij ( |2005[) fitting 
results are consistent with this. Although our 3 Gyr, 0.4 Z© popu- 
lation only accounts for 0.12 % of the stellar mass contributing to 
the spectrum, it do minates the flux at w avelengths < 4500 A. 

More recently, Buson et al.l ( f2004h . by matching the HST/FOS 
UV spectra of NGC 5018 with that of the well-studied compact 
dwarf elliptical M32, claimed that NGC 5018 is a galaxy whose 
central regions are dominated by a 3 Gyr stellar population, and 
that therefore, this galaxy is the child of a major gas-rich merger. 
They ask the question: 'is there any evidence of an older stellar 
population in NGC 5018, or were literally most of the stars we see 
formed in the merger event?'. Using our long-baseline fitting tech- 
nique, we can of course answer that question. We do indeed find a 



mature stellar population, which dominates the total stellar mass, 
but, as the young population dominates the flux in the UV, this can- 
not be detected without using the full spectral range. Conversely, of 
course, we cannot detect the younger population without including 
the UV wavelengths. 

Although the one- and two-component fitting algorithms do 
not consider any dust extinction of the spectrum, it is clear from 
the MOPED results that this is not a large effect in either galaxy. 
The quality of the fit over an extended baseline endorses the re- 
covered dust parameters, as an error in the dust estimation would 
result in a gradual offset of the fitting continuum over the length of 
the spectrum. The one- and two-component fits suffer less in this 
respect, as the relative normalisations of the sections are allowed to 
vary. 

From the detailed fitting, we do not recover formally accept- 
able v alues of . Recent work (e.g. lThomas, Maraston & Bended 
l2003b has found super-solar a-element abundances present in 
early-type galaxies. With non-solar abundances calculated for in- 
dividual features, they are able to recover better and more consis- 
tent fits to absorption-line features than for the solar abundances. 
Non-solar a-elements probably account for individual regions of 
poor fit to the data in this work, but, as long-baseline, non-solar 
Q-element stellar population models do not as yet exist, we cannot 
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object 


Te 1 kpc 


He (K / mag) 


(TO / kms ^ 


NGC 3605 


1.9 


17.8 


103 


NGC5018 


4.3 


16.6 


223 


giant core 


7.54 ± 1.70 


17.36 ± 0.68 


238 ±1 


power-law 


2.29 ± 1.75 


16.75 ± 1.00 


157 ±2 



oo 



Table 4. Observed properties of NGC 3605 and NGC 5018: tro is the central 
velocity dispersion; /i^, the effective surface brightness in the K band, is 
calculated from the mean surface brightness within the effective half-light 
radius (re) in the B band, using (B — K) = 3.9 (Genzel et al.. 2001). re 
for NGC 5018 is from Scorza et al. (1998), re for NGC 3605 is from Faber 
et al. (1997), and /Je and ctq for both galaxies is from Faber et al. (1989). 
Also shown are the mean properties of the galaxy samples plotted in Figure 
1101 The limits shown are the standard deviation of the sample. All data are 
taken from Bender et al. (1992), Faber et al. (1997) and Faber et al. (1989). 



test whether the quality of the fit will be significantly improved by 
their use. 

Even though, for both galaxies, the vast majority (> 99%) of 
the stars are reproduced by only two stellar populations, the unre- 
stricted star-formation history permitted by the MOPED algorithm, 
together with its inclusion of a dust parameter, has enabled a bet- 
ter quality of fit to the two low-redshift elliptical galaxy spectra 
than that which is possible with even the best-fitting two-population 
model, i.e. for NGC 3605, x2= 1-8 (compared with 2.1 for the two- 
component model), and for NGC 5018, xi= 1.4 (also 2.1 for the 
two-component model). 

Despite all of the details discussed above, it is encouraging 
that, the detailed spectral fits, as well as the MOPED fits, are in 
good agreement regarding the ages of the populations in these 
galaxies. 

5.2 Mergers and the fundamental plane 

Figure [To] shows projections in the fundamental plane of the ob- 
served properties ofrTOCJ6^^ and NGC 
5018 jpaber et al.iri989l : IScorza et al.lll998h . together with a popu- 
lation of giant boxy ellipticals, and intermediate mass, disky ellip- 
ticals jBender e t al~ 1992; Faber etal. 1997) for comparison. The 
two galaxy populations occupy different positions in the funda- 
mental plane. Giant ellipticals, with core central surface bright- 
ness distributions, tend to be more massive, more luminous and 
have higher central velocity dispersions than the intermediate-mass 
ellipticals with power-law central surface brightness distributions 
('cuspy' galaxies). Mean values of the effective surface brightness, 
/Xe, effective radius, re, and central velocity dispersion, a, for the 
two galaxy populations are listed in Table |4l together with the cor- 
responding values for NGC 3605 and NGC 5018^ 

The HSTAVF PC2 observa t ions o f both iFaber et al.1 h997t) 
and, more recently, iLauer et ^ j2005l) , reveal that NGC 3605 is 
a cuspy galaxy, and its position on the fundamental plane places 
it unequivocally in the cuspy galaxy population. The spectroscopic 
determination of the ages of the constituent populations of NGC 
3605 is consistent with the theory that cuspy galaxies are the prod- 
uct of a recent gas-rich galaxy-galaxy merger. The bulk of the stars 
are old (> 10 Gyr), but there is also a significant young (1 Gyr) 
stellar population, with a higher metallicity. The higher metallicity 
is expected to result from the burst of star formation triggered in a 
recent major gas-rich merger, where the new stars are formed from 
gas which was metal-enriched when the older stars were formed.. 

MOPED cannot at present fit multiple spectra simultaneously. 
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Figure 10. Projections of the fundamental plane. Green squares; giant / 
core ellipticals (from Bender et al., 1992, and Faber et al., 1997); black 
circles; intermediate mass / power-law ellipticals (from Bender et al., 1992, 
and Faber et al., 1997); red star; NGC 3605 (from Faber et al., 1989 and 
Faber et al., 1997); blue triangle; NGC 5018 (from Faber et al., 1989 and 
Scorza et al., 1998). All parameters are calculated for a cosmology with 
Ho = 70 kms~^Mpc~^ and f2m = 0.3. See section 5.3 for discussion. 



In this case, therefore, adjoining spectral segments were normalised 
to unity in the overlap regions, and spliced together. 

NGC 5018, a shell galaxy, is at first sight more ambiguous. 
Its effective radius, central velocity dispersion and relatively low 
fie suggest that it belongs with the giant core galaxy population, 
as one might expect from a brightest group elliptical. At the same 
time, NGC 5018 appears to contain a fairly substantial young pop- 
ulation (12% by mass). One might expect this young population 
to have been triggered by a merger event. However, NGC 5018 is 
a near companion to the spiral galaxy NGC 5022, and there is evi- 
dence of an on-going interaction between thes e two galaxies: an op- 
tical b ridge has been detected bet ween them jSchweizer & Seitzej 
Il992h , together with an HI bridge jKim etalj 19881) . Kim et al. con- 
clude that the HI interaction between NGC 5018 and NGC 5022 
is current, but ,Malin & Hadlev ( 1997) infer from the deeply em- 
bedded dust lane in NGC 5018 that interactions between the two 
galaxies have occurred more than once in their history. Simulations 
indicate that shell structure resulting from a weak interaction be - 
tween galaxies may survive for a few Gyr ( lReduzzietid]|l996h . 
Therefore, although extra blue light from the shells and bridge 
mimics a post-merger galaxy, it is more likely that NGC 5018 
consists of a predominantly mature central population, which has 
evolved passively, together with younger shells which are a prod- 
uct of mass transfer or accretion. In addition, the lower metallicity 
of the younger population (0.4 Z©) compared with the older pop- 
ulation (2.5 Zq) suggests that the more-recent stars were formed 
from gas which is either not metal-enriched, or has been diluted by 
un-enriched gas, which is consistent with star formation triggered 
by the infall of HI along the bridge between NGC 5018 and NGC 
5022. The position of NGC 5018 in the fundamental plane favours 
this scenario over that of a major merger between two early-type 
galaxies. 
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6 CONCLUSIONS 

We have undertaken a detailed exploration of the star-formation 
histories of the elliptical galaxies NGC 3605 and NGC 5018. By 
assembUng and analysing the spectral energy distributions of these 
galaxies spanning the wavelength range 2500 - 8000 A we aimed 
to i) determine the added value of the ultraviolet extension for 
breaking the age-metallicity degeneracy, and disentangling star- 
formation history, and ii) properly constrain and compare the star- 
formation histories of these two elUpticals, which occupy very dif- 
ferent regions of the fundamental plane. Multi-component stellar 
population fitting allows us to uncover a far more complete picture 
of the star-formation history of early-type galaxies than fitting to 
absorption line features. 

We find that i) optical spectra with A > 3500 A may not con- 
tain sufficient information to robustly uncover all the stellar popula- 
tions present in individual galaxies, even in such relatively passive 
objects as elliptical galaxies, ii) the addition of the ultraviolet data 
approaching A = 2500 A holds the key to uncovering the major 
epochs of star formation for both galaxies, and establishing well- 
constrained star-formation histories, from which we can infer a for- 
mation and evolution history which is consistent with their photo- 
metric properties. Our long-baseline fitting technique is remarkably 
robust to choice of stellar population model and wavelength reso- 
lution iii) despite the superficial similarity of their spectra, the two 
galaxies have very different star-formation histories - the smaller, 
cuspy elliptical NGC 3605 contains a high-metallicity population 
of age 1 Gyr, and has a position on the fundamental plane typical 
of the product of a low-redshift gas-rich merger, while the giant el- 
liptical NGC 5018, with a sub-solar secondary population, appears 
to have gained its more recent stars via mass transfer / accretion of 
gas from its spiral companion, iv) despite these differences in de- 
tailed history, more than 85% of the stellar mass of both galaxies is 
associated with an old (9-12 Gyr) stellar population of near-solar 
metallicity. 

For a universe with fim = 0.3, Oa = 0.7 and Ho = 70 
kms~^Mpc~^, it thus appears that the redshift at which NGC 3605 
was assembled to its present state is Zmerg = 0.076, if the for- 
mation of the stars in the 1 Gyr population is assumed to have been 
triggered by this event. This is in agreement with hierarchical struc- 
ture formation scenarios in which disk-disk merging events take 
place predominantly at z < 1. However the vast majority of its stars 
appear to have formed at Zf > 3. NGC 5018 also formed the vast 
majority of its stars at high redshift (z/ = 2 — 5), but then appears 
to have gained a shell (or shells) of younger stars between z ~ 0.3 
and the present-day, possibly through interaction with NGC 5022. 

In conclusion, the spectroscopic disentanglement of the com- 
ponent populations of elliptical galaxies is possible, and is a use- 
ful tool in the investigation of elliptical galaxy formation and evo- 
lution. We are able to determine the formation scenarios of these 
galaxies, and hence their general position in the fundamental plane 
using this method. The availability of full SEDs with non-solar a- 
element abundances would help to further improve the goodness- 
of-fit achieved with this method. 



ACKNOWLEDGEMENTS 

This work was based in part on observations with the NASA/ESA 
Hubble Space Telescope, obtained at the Space Telescope Sci- 
ence Institute, which is operated by the Association of Universi- 
ties for Research in Astronomy, Inc. imder NASA contract No. 



NAS5-26555. Louisa Nolan acknowledges the support of PPARC, 
through the award of a PDRA. Raul Jimenez is supported by 
NSF grants AST-0408698 and PIRE-0507768, and NASA grant 
NNG05GG01G. Ben Panter is supported by the Alexander von 
Humboldt Foundation, the Federal Ministry of Education and Re- 
search, and the Programme for Investment in the Future (ZIP) of 
the German Government. 



REFERENCES 

Archibald, E. N., Dunlop, J. S., Hughes, D. H., Rawlings, S., 

Eales, S. A., Ivison, R. J. 2001, MNRAS, 323, 417 
Baugh CM., Cole S., Frenk C.S., 1996, MNRAS, 283, 1361 
Baugh CM., Cole S., Frenk C.S., Lacey C.G., 1998, ApJ, 498, 
504 

Begelman, M. C, Blandford, R. D., Rees, M. J. 1980, Nature, 287, 
307 

Bender R., Surma P., Dobereiner S., Mollenhoff C, Madejsky R., 

1989, A&A, 217, 35 
Bender R., Burstein D., Faber S.M., 1992, ApJ, 399, 462 
Beuing J., Dobereiner S., Bohringer H., Bender R., 1999, MN- 
RAS, 302, 209 
Bica E., AUoin D., 1987, A&A, 186, 49 
Binney J., 2003, MNRAS, 347, 1093 
Bruzual, G., Chariot, S. 2003, MNRAS, 344, 1000 
Buson, L. M., Bertola, F., Bressan, A., Burstein, D., Cappellari, 

M. 2004, A&A, 423, 965 
Caldwell N., Rose J.A., Concannon K.D., 2003, ApJ, 125, 2891 
Cimatti, A., et al. 2004, Nature, 430, 184 

De Lucia, G., Springel, V., White, S. D. M., Croton, D., & Kauff- 

mann, G. 2006, MNRAS, 366, 499 
Dorman B., O'Connell R. W., Rood R. T, 2003, ApJ, 591, 878 
Dunlop J.S., et al, 1996, Nature, 381, 581 
Dunlop J.S., McLure R.J., Kukula M.J., Baum S., O'Dea CP., 

Hughes D.H., 2003, MNRAS, 340, 1095 
Eisenstein D.J., et al., 2003, ApJ, 585, 694 
Faber S.M., et al., 1997, AJ, 1 14, 1771 

Faber, S. M., Wegner, G., Burstein, D., Davies, R. L., Dressier, A., 

Lynden-Bell, D., Terlevich, R. J. 1989, APIS, 69, 763 
Faber S.M., et al., 1997, AJ, 114, 1771 

FaneUi M.N., O'Connell R.W., Burstein D., Wu C.-C, 1990, ApJ, 
364, 272 

Ferrarese, L., Merritt, D. 2000, ApJL, 539, L9 
Fioc M., Rocca-Volmerange B., 1997, A&A, 326, 950 
Gallagher J. S., Hunter D. A., Bushouse H., 1989, AJ, 97, 700 
Genzel R., et al., 2001, ApJ, 563, 527 

Gordon K. D., Clayton G. C, Misselt K. A., Landolt A. U., Wolff 

M. J., 2003, ApL 594, 279 
Heavens, A.F., Panter, B.D., Jimenez, R., Dunlop, J.S., 2004, Na- 
ture, 428, 625 

Jansen, R. A., Fabricant, D., Franx, M., Caldwell, N. 2000, ApJS, 
126, 331 

Heavens, A. F, Jimenez, R., Lahav, O. 2000, MNRAS, 317, 965 

Kauffman G., 1996, MNRAS, 281, 487 

Kauffmann, G., & Chariot, S. 1998, MNRAS, 294, 705 

Kauffman G., Chariot, S. 1998, MNRAS, 297, L23 

Kauffman G., et al., 2003, MNRAS, 341, 54 

Khochfar, S., & Burkert, A. 2005, MNRAS, 359, 1379 

Kim, D.-W., Jura, M., Guhathakurta, P., Knapp, G. R., & van 

Gorkom, J. H. 1988, ApJ, 330, 684 
Kormendy J., Bender R., 1996, ApJL, 464, LI 19 



Lauer T.R., et al., 1995, AJ, 1 10, 2622 

Lauer, T. R., et al. 2005, AJ, 129, 2138 

Khochfar, S., & Burkert, A. 2005, MNRAS, 359, 1379 

Lotz J.M., Ferguson H.C., Bohlin R.C., 2000, ApJ, 532, 830 

Malin, D., & Hadley, B. 1997, Publications of the Astronomical 

Society of Australia, 14, 52 
Maraston, C, 2005, MNRAS, 362, 799 
McCarthy, R J., et al. 2004, ApJL, 614, L9 
Panter, B., Heavens, A. R, & Jimenez, R. 2003, MNRAS, 343, 

1145 

Panter, B., Jimenez, R., Heavens, A. F & Chariot, S., 2006, sub- 
mitted to MNRAS, astro-ph/0608531 
Ponder, J. M., et al. 1998, AJ, 116, 2297 

Ravindranath S., Ho L.C., Filippenko A.V., 2002, ApJ, 566, 801 
Reduzzi, L., Longhetti, M., & Rampazzo, R. 1996, MNRAS, 282, 
149 

Reichardt, C, Jimenez, R., Heavens, A. F. 2001, MNRAS, 327, 

849 

Renzini, A., 2006|astro-ph/06 03479 to appear on Annual Review 

of Astronomy & Astrophysics, Vol. 44 
Scorza, C, Bender, R., Winkelmann, C, Capaccioli, M., & Mac- 

chetto, D. F 1998, A&AS, 131, 265 
Scott, S. E., et al. 2002, MNRAS, 331, 817 
Schweizer, F, Seitzer, P 1992, AJ, 104, 1039 
Silk, J., Rees, M. J. 1998, A&A, 331, LI 
Stevens, J. A., et al. 2003, Nature, 425, 264 
Terlevich, A. I., Forbes, D. A. 2002, MNRAS, 330, 547 
Thomas, D., Maraston, C, Bender, R. 2003, MNRAS, 339, 897 
Trager S.C., 1999, In: Photometric Redshifts & High-z galaxies, 

ASP Conf Ser vol 191, p. 195 
Trager S.C., et al., 2001, In: Astrophysical Ages and Time Scales, 

(a stro-ph/0 1 04404) 
Worthey G., 1994, ApJS, 95, 107 
Worthey, G., Ottaviani, D. L. 1997, ApJS, 111, 377 



SFH of ellipticals across the fundamental plane 



